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We calculate analytically a back reaction of the stationary spherical accretion flow near the event 
horizon and near the inner Cauchy horizon of the charged black hole. It is shown that corresponding 
back reaction corrections to the black hole metric depend only on the fluid accretion rate and 
diverge in the case of extremely charged black hole. In result, the test fluid approximation for 
stationary accretion is violated for extreme black holes. This behavior of the accreting black hole is 
in accordance with the third law of black hole thermodynamics, forbidding the practical attainability 
of the extreme state. 
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I. INTRODUCTION 

The matter and fields, inflowing into a rotating or 
charged black hole, generate the drastic singular phenom- 
ena at the inner Cauchy horizon: the infinite blue-shift 
and mass inflation (see e.g. |ll-[l9|). We describe here a 
back reaction of the stationary spherical accretion flow 
and reveal the divergence of back reaction perturbations 
at both the event horizon and the inner Cauchy horizon 
of the nearly extreme charged black hole. This singu- 
lar behavior of back reaction perturbations implies the 
violation of test fluid approximation for the stationary 
spherical accretion onto extreme black holes. 

The stationary accretion onto black hole is a solution 
with an independent of time t and radius r influx of test 
fluid (or some field) with the energy density poo and 
pressure Poo far from the black hole. We describe be- 
low the corrections to the background spherically sym- 
metric metric of the black hole due to the back reaction 
of stationary accreting perfect fiuid. The problem of a 
stationary accretion onto black hole is self-consistent, if 
(i) the accreting fiuid is lightweight, and (ii) the rate of 
black hole mass growth is slow (the stationary limit). 
Two small parameters are needed to satisfy these condi- 
tions. The first small parameter is a mass ratio of the 
inflowing gas in the volume with the black hole gravita- 
tional radius, rrigas ~ poo'm^, to the black hole mass m: 
PoaiTi^ /m — PooTT^^ ^ 1- With this small parameter, the 
test fluid approximation in the background metric is valid 
inside the region r <^ i?max — 'm{pryom?)~^'^ around the 
black hole. The second small parameter is a slowness 
of the black hole mass changing with respect to the hy- 
drodynamical time, m/m <C Cs/m, where Cs is the sound 
velocity in the accreting gas. For example, the stationary 
spherical accretion rate of the perfect fluid is [20, l2l[ : 



TO = 47rv4(poo +Poo)"l^ 



(1) 



where A is a numerical constant, depending on the equa- 
tion of state p = p{p)- This constant A is of the order 
of unity for relativistic fluids with p ^ p {i. e, for flu- 
ids with a relativistic sound velocity Cs = 0{1)). Two 



the relativistic perfect fluids. The back reaction of fluid 
on the background black hole metric in the considered 
stationary accretion limit may be found by perturbation 
method due to existence of these two small parameters. 
Note that in the stationary accretion limit the apparent 
horizon will coincide with the event horizon of the ac- 
creting black hole. Quite the different approaches were 
used in |23 - [25j for analysis of the back reaction on the 
Reissner-Nordstrom black hole. 

Below we describe the back reaction of the stationary 
accretion with a linear accuracy with respect to the small 
parameter to = const <C 1. We find analytically the 
corresponding back reaction corrections to the black hole 
metric metric near the event horizon and near the inner 
Cauchy horizon. 



II. EINSTEIN EQUATIONS 

A spherically symmetric metric may be written in the 
general form with two arbitrary functions i26i : 



ds^ = e'^^'^^Ut^ ~ e^(*''-)dr2 



'-{de^ 



ld02). (2) 



For application to the back reaction of the accreting 
matter on the charged black hole, we define two metric 
functions, fo{t,r) and fi{t,r), and also two correspond- 
ing "mass functions" (or mass parameters) toq {r, t) and 
TOi(r,t): 



s''^*'''^ - /o(t,r) = l- 



-\{t,r 



' = h{t,r) = l 



2mo(t,r) 

r 
2mi{t, r) 



(3) 
(4) 



small parameters become identical, to 



< 1, for 



where e is the black hole electric charge. In the case of 
the Reissner-Nordstrom metric (i. e., in the absence of 
accreting fiuid), both mass functions coincide with the 
black hole mass, ttiq = to,i = to,, and, respectively, /q = 

A spherically symmetric gravitational field in the gen- 
eral case is defined by the four Einstein equations. Three 
of them are the differential equations of the first-order. 



and the fourth one is of the second-order. These equa- 
tions for metric © have the following form [26| : 



SttTJ = -e 

SnT^ = 

St:TI = 

8ttT^ = 
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v" -F {v' - A') TT + 



(5) 
(6) 

(7) 

(8) 

(9) 



The corresponding components of the energy-momentum 
tensor for a perfect fluid are 
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(11) 
(12) 
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where u = dr/ds is the radial component of the fluid 4- 
velocity, and, respectively, p and p is an energy density 
and pressure of fluid in the comoving frame. Below it is 
supposed an arbitrary equation of state p = p{p), relating 
the fluid pressure and energy density. The Bianchi iden- 
tity holds true for Einstein equations, and so only three 
equations from four in pU)) -([T ^ are independent. We 
choose ([U, dU and ([7]) for these independent equations. 

Zero approximation in our approach corresponds to the 
stationary spherically symmetric inflow of the test fluid 
in the background Reissner-Nordstrom metric. The cor- 
responding solution [20, [2l|, [231 defines the conserved ra- 
dial flux of energy rh and the radial dependance for the 
4- velocity component u = dr/ds = ""(r), for energy den- 
sity p = p(r) and for pressure p — p(p) — p(r). Respec- 
tively, this solution fixes all components of the energy- 
momentum tensor T*'^. For self-consistency of the accre- 
tion problem in the background metric, the radial infiux 
of energy must be small, i. e. m ^ 1. 

In the first approximation we will take into account 
the linear contributions with respect to m ^ 1 to the 
energy momentum-tensor in Einstein equations. As a re- 
sult, we find the deviation of metric from the background 
one, i. e. the back reaction with a linear accuracy with 
respect to the small parameter m. It will be shown that 
the back reaction corrections near the event horizon and 
inner Cauchy horizon of the black hole depend only on 
the accretion rate m and do not depend on the equa- 
tion of state p = p{p) of the accreting fluid. Formally, to 



find a back reaction modification of the black hole metric 
near the horizons, it is needed to consider the space-time 
region, where both /o ^ 1 and /i <C 1. 



III. BACK REACTION IN THE 
SCHWARZSCHILD METRIC 

As the first step, we find the back reaction of accretion 
near the event horizon of the Schwarzschild black hole. 
The first Einstein equation ([SJ defines the conserved ra- 
dial flux of anergy, i. e. the matter accretion rate. 



= -A^r\p+p)uJ-^{h 



const. 



(14) 



The value of this flux ([T]) is defined (in the zero ap- 
proximation) from the stationary solution of the test 
fiuid accretion in the background Schwarzschild metric 
[13 HH, m m. *** The radius of the modified event 
horizon, r = r+, in the considered stationary accretion 
limit is defined by the condition /o(r_|_,t) = 0. The Ein- 
stein equations are hyperbolic in general, and so we have 
at the horizon the second condition, fi{r^,t) — 0. *** 

The first Einstein equation ([5|) in the linear approxi- 
mation with respect to to <C 1 has a very simple form. 



drrii 



(15) 



This equation (in the considered stationary accretion 
limit with to ^ 1) alludes the "factorization" of the 
mass functions: mQ{r,t) — m{t) pQ(r) and, respectively, 
mi{r,t) = m{t) pi{r). The dimensionless mass functions 
Po{f) and pi{r) are defined so, that fJ,o{r) — Pi{r) = 1 at 
TO = 0. After substitution of this "factorization" ansatz 
in (fTSl). we obtain 



dm{t) 
dt 



Pi{r) 



(16) 



The r.h.s. of this equation is already linear with respect 
to TO ^ 1. So, for the dimensionless function pi{r) in this 
equation it is needed to use a zero approximation with 
respect to to, i. e. to put pi{r) = 1. In result, the partial 
differential equation (llSp for the function toi (r, t) reduces 
to the ordinary differential equation (|16p for m{t). The 
corresponding solution of this reduced equation is 



m{t) = 77i(0) + / m{t')dt'. 



(17) 



Here to.(0) is a black hole mass at some initial moment 
t = and 7n(t) is a current value of black hole mass. 
The black hole mass m{t) is slowly growing due to a slow 
accretion rate to <C 1. This is the only one equation 
in our analysis, where the temporal dependance of the 
accretion rate m{t) must be taken into account. In all 
other equations the accretion rate m{t) is only a small 
constant parameter. 



At this step we may find the function /xi(r) with the 
help of the second Einstein equation ([5]). After factor- 
ization of the mass functions ttiqI and mi, it is useful 
in the following to use the dimensionless radial variable 
X — r/m{t) with m{t) from p7|) . The second Einstein 
equation ([6]) is written now in the form 



ax 



P+{p + p)-r 
n 



(18) 



A combined solution of equations p^ and ((T5| with the 
using of accretion solution for fluid with equation of state 
p = p{p), defines the requested function /ii. Near the 
event horizon, where /o ^ 1 and /i <C 1, from (1141 and 
P^ we obtain 



dx 



Attx^ {p + p)v? 



2™ 



2^1 (x) 



(19) 



By using a new variable (5(x) = x — 2/ii(a;) ^ 1, we 
find the approximate solution of equation (J19l) near the 
event horizon x+, defined by the condition 5{x+) — x+ — 
2^l^{x+) = Q■. 



HiiS) sa /x+ + 2m log 



Am 



where 



/i+ ~ Pi{x+) w 1 + 2m log |m|. 



(20) 



(21) 



is the value of mass function /ii(a;) at the event horizon 
radius x+, modified by the back reaction. Solution ((20)) 
for the inverse function x = x{fii) is written in form: 



x{fii) sa 2fi++4m 1 + 



2m 



exp ■ 



Hi - fi. 



2m 



, (22) 



where x+ = x{fi+) = 2/x+. Solution (|20l) for the 
Schwarzschild metric, modified by the back reaction of 
accreted matter is shown in Fig. 1. It is important to 
note that this solution is valid only in the narrow region 
8m = |(5min| < |(5| ^ 1 around the modified event hori- 
zon x+. The used linear approximation with respect to 
m would be insufficient for calculation of the correspond- 
ing mass function at \S\ < |(5min| (inside the filled box in 
Fig. 1). 

During integration of equation (J18p we retained only 
the leading term with /i ^ 1 in the denominator, provid- 
ing the major contribution to the solution ^ mlog|m|, 
and neglected all contributions to the solution of the or- 
der of m. Formally, we supposed that distribution of fluid 
around the black hole is a finite sphere of radius Xq , sat- 
isfying the condition 1 ^ Xq ^ Xmax = (PooW^)"^/'^. 
We neglect the contribution to solution of equation (fT5)) 
from the term 4ttx^p, related with mass of the accreting 
gas inside the sphere of radius Xq, and also the contri- 
bution from the leading term at the upper integration 
limit ^ m/Xo, which is the gravitational "mass defect". 
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FIG. 1: The mass function /ii(x) from equation (|20|l near the 
event horizon of black hole x+, modified by back reaction, 
in the linear approximation with respect to rh <S^ 1. Inside 
the filled box the used linear approximation is insufficient for 
determination of the mass function pi{x). 



Finally, we neglect in the integration constant the contri- 
bution from the boundary condition, which is also ~ m. 

The third Einstein equation ([7]) for the considered per- 
fect fluid has the form 

„,2 



/i ,^1 
Jo X 



Ml 
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Mo-p 
Jo 



-47ra; 



(p + p)-r+p 

/i 



Near the event horizon this equation is written as 

-,,2 



Mo 



-47ra; (p + p) 



h 



-Ml- 



(23) 



(24) 



Here it is taken into account that a ratio /o//i near the 
horizon equals to its background value, (/o//i)+ = 1, at 
the linear approximation with respect to m ^ 1. So- 
lution of equation (j24l) near the event horizon, where 
6 ^ X - 2fii{x) < 1, with the help of ^^ and ([24]) 
is 



/io((5) « p+ 



2771 log |1-— I 
4m 



(25) 



By comparing (PD)) and (pS)) . we see that near the horizon, 
at I (5 1 ^ 1, it is satisfied the condition 



Mo('5) + Mi(^) ~ 2^^ 



(26) 



To calculate the corresponding back reaction in the case 
of the Reissner-Nordstrom black hole ([2]) - (|1]) we define 
the extreme parameter of the black hole e = y^l — e^ /ra^ 
and a new variable 5± = x— [p-i{x) ± •\/mi(^)^ — 1 + e^]. 




FIG. 2: The shallow funnel is the embedding diagram for the 
Shwarzschild BH (Flamm's paraboloid) and the deep funnel 
is the corresponding embedding diagram for the near extreme 
Reissner-Nordsrom BH. The height of the funnel tends to in- 
finity at e — >■ according to (pT|). 



IV. BACK REACTION IN THE 
REISSNER-NORDSTROM METRIC 

Quite similar to the Schwarzschild case, near the event 
horizon a;+ and the inner Cauchy horizon x^ of the 
Reissner-Nordstrom black hole we find: 



Hli5±) r 


(\±ef . , 
:s l±^_-^mlog 


Mo('5±) ' 


- 2y^± -Mil^i), 


where 
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-l±(^±^)%.lo. 
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m — S± 



2e 



2e2 



,(27) 
(28) 

(29) 



This solution is vahd only at [(1± e)^/e'^]m < |(5±| <C 1 
and (5_ > 0. The value of fii{x) at the modified event 
horizon x = x+ and at the inner Cauchy horizon x = x^ 
corresponds formally to S± — in (|27p . The correspond- 
ing radii of horizons are 



x± 



Ml 



ix±) ± y/fi,ix±y-il-e^) 



(30) 



(lie) 



-K^ 



mlog 



(l±e)3 



We uphold here the major back reaction term ex m log \m\ 
and neglect the much smaller term ex to. 



if e — >■ 0. This behavior is in agreement with the cos- 
mic censorship conjecture [29[ and with the third law 
of black hole thermodynamics [3^ : the extreme state is 
unattainable in the finite processes or, in other words, 
it is impossible in practice to transform the black hole 
into the naked singularity. A similar conclusion was de- 
rived recently in [33] by demonstration that the Reissner- 
Nordstrom black hole can never be overcharged to the 
naked singularity via the absorption of a charged parti- 
cle. 

Note, that the test fluid approximation is violated for 
the stationary accretion of ultra-hard fluid with p = p 
at the event horizon of the extreme Kerr-Newman black 
hole even without the back reaction corrections [21, 3l|, 
|32| . Violation of the test fluid approximation for the 
accretion onto extreme black holes is also in accordance 
with the absence of solutions for stationary accretion onto 
the naked singularities [2l|, |32| . To resolve the problem 
of back reaction of accreting matter on the extreme black 
hole it is requested to find a solution of Einstein equations 
beyond the perturbation level. 

A physical reason for the divergence of back reaction 
corrections in the extreme case is in the infinite stretching 
of a local distance l{x) to the event horizon r_|_ (and, 
analogously, to the Cauchy horizon r_) at e — >■ 0: 
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See in Fig. [5] the corresponding embedding diagram for 
the Reissner-Nordsrom BH, constructed from the relation 
dz'^ + dx^ = f'^dx^ with / from ^ and (g]). The space 
near horizons is loaded with a finite energy density p{x) of 
the accreting fluid along the local distance l{x), providing 
contribution to the mass functions m(){x) and mi{x). 

Due to the infinite stretching of l{x) in the extreme 
case e = 0, the total mass of accreting gas and the cor- 
responding mass functions ij,o{x) and Hi{x) in ((27|) and 
(1^51) are diverging near the extreme black hole. Quite a 
similar infinite stretching of l{x) is also a characteristic 
feature of the extreme Kerr black hole [Sj. Therefore, 
the described violation of the test fiuid approximation 
for the stationary accretion might be crucial both for the 
extremely charged and extremely rotating black holes. 



CONCLUSIONS 



From ^7} - (1201) it follows directly, that the test fluid 
approximation is violated due to the back reaction of the 
accreted fluid in the extreme black hole limit e — > 0. 
Namely, the corresponding corrections to the radius of 
the black hole event horizon and the inner Cauchy hori- 
zon diverge at any arbitrarily small accretion rate to. 
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